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Abstract 



We propose a Hamiltonian model that describes the interaction between 

a vortex line in superfluid 4 He and the gas of elementary excitations. An 

equation of irreversible motion for the density operator of the vortex, regarded 

as a macroscopic quantum particle with a finite mass, is derived in the frame 

of Generalized Master Equations. This enables us to cast the effect of the 

coupling as a drag force with one reactive and one dissipative component, in 

agreement with the assumption of the phenomenological theories of vortex 

mutual friction in the two fluid model. 
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I. INTRODUCTION 



Since the discovery of quantized vortices in liquid helium II, it has been recognized that 
they might provide a mechanism for the coupling of the superffuid to the normal fluid. In 
the two fluid model, this coupling is represented by a mutual drag force with one dissipative 
and one conservative component, whose respective strengths can be measured investigating 
the attenuation of second sound at various temperatures |j]-BJ. Models for the friction 
coefficients which successfully fit the data up to 2.1 K have been presented in Refs. [|]||. The 
phenomenon of vortex mutual friction has been observed as well in rotating superffuid 3 He-B 
H; this fact brings support to the conjecture that such a mechanism is indeed a relevant 
source of dissipative processes in highly degenerate quantum fluids. It is also worthwhile 
to remind that in the last years, vortices have been seen to play a role in phase transitions 
taking place either in underpressurized 4 He or in supersaturated solutions of 3 He in 4 He 
i- 

The interaction between the velocity field and the density fluctuations of the superfluid, 
which at nonvanishing temperatures are embodied in the normal fluid, can be accounted for 
within a Lagrangian description fl,|T0|| . However, microscopic descriptions of the interaction 
between the superfluid motion, especially when topological singularities are concerned, and 
collective excitations, are not conclusive The investigation of quantum tunneling of 

vortex lines in superconductors and superfluids fll~HT3| has now improved our comprehension 



of, for example, the role of the inertial mass of the vortex (see also Refs. [|14|-|16[| and cited 
therein) and the influence of either pinning or dissipation on the tunneling rates. Indeed, 
the value of the vortex inertia is a fundamental parameter in any theoretical description 
of vortex dynamics and it remains being a controversial issue |16|J17| . Different starting 
points assign to this inertia figures ranging from zero to infinity, the latter arising from 
a logarithmic divergence with the system size due to the renormalization effect induced 
by the condensate motion. On the other hand, the vortex mass is known to be finite in 
superconductors [pjfl . Although the phenomenological approaches completely disregard 
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inertial effects, the mass enters the description of the dynamics of a free vortex, known to 
be cyclotron-like [Q,|l(],|l7]] , through a frequency parameter Q, which would be a measurable 
quantity if this time dependent regime were experimentally visualized. In particular, it has 
been recently shown [K| that the cyclotron motion is a natural solution of the nonlinear 
Schrodinger equation applied to a vortex. Dynamical |TIJ and thermodynamical |HJ methods 
have been proposed to measure the vortex inertial coefficient; our present purpose is neither 
to participate in the existing polemics, nor to propose a new model for the calculation of 
the vortex mass, but rather to assume that it is a numerical parameter and proceed along 
similar lines as those invoked in the well established cyclotron motion already discussed in 
textbooks 0. 

The aim of the present work is to propose a Hamiltonian model for the coupling between 
a rectilinear vortex immersed into the excitations of the superfluid, as will be discussed in 
Sec. 2. Due to the translational symmetry of the problem along an axis parallel to the 
vortex, the problem is a twodimensional one, i. e., we consider a point vortex on a plane. 
We shall show in Sec. 3 that, if one considers the vortex as a quantum particle undergoing 



Brownian motion EH] in a heat reservoir, it is possible to establish the irreversible time 



evolution of its density operator within the Generalized Master Equation (GME) approach 



21|j22|| . In this way, in Sec. 4 we are able to derive dissipative equations of motion for the 



canonical position-momentum variables of the vortex and for its velocity. These variables 
can be seen to evolve under the combined effect of the usual hydrodynamical lift on a rotating 
cylinder, plus a drag force. If the coupling is linear in the excitation operators, the drag 
coefficients are governed by the dynamical susceptibility of the liquid. The consequences of 
the equations of motion thus obtained, the asymptotic velocity of the vortex and the relation 
of the current description to the phenomeno logical model are discussed in Sec. 5, where the 
perspectives of the present approach are also outlined. 
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II. THE HAMILTONIAN MODEL 



Let us first summarize the description of the free motion. The Hamiltonian for a cylin- 
drical vortex parallel to the z— axis in liquid helium at zero temperature is in charge of 
providing the Magnus force 0. It reads 

Hv = Wi [p " gA(r)]2 + MQv *y 

where 

A(r) = ^(y,-x) (2.2) 

is the vector potential whose curl yields the vortex-velocity-dependent part of the Magnus 
force and the potential term MQv s y gives the superfluid-velocity-dependent part of this 
force. Here M is the dynamical mass of the vortex, p s denotes the number density of the 
superfluid, v s its velocity along the assumed to be uniform, h is Planck's constant, 

/ the system length along the z— axis and 

n = £^i. (2.3) 

The quantity q — ±1 is the sign of the vorticity according to the right handed convention. 
Furthermore, at zero temperature, p s coincides with the total density per unit mass p/m, 
being m the mass of a helium atom. 

At this point it is convenient to remember the existing theoretical uncertainty regarding 
the vortex mass parameter M that should appear in dynamical calculations [T^|l^JT^] and 



keep in mind that in the phenomenological approaches || the dynamical regime of the vortex 



is that in which the Magnus force balances the drag plus any applied force [23|. As stated in 



the Introduction, our viewpoint here is identical to that of former authors @,|II],|Tl],|n| who 
assume a finite figure for the vortex inertia and consider the cyclotron - like motion of a free 
vortex as their starting point, with the frequency Q as the leading parameter. Since it will 
be shown in Sec. 4 that the dissipative motion is easily described in terms of the complex 
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position variable z = x + iy and the velocity dz/dt, we here write the complex Hamilton 
equation that stems from (2.1) 

d 2 z ( dz \ * v 

^ =lQ [*-"•) {2A) 

with the complex Magnus force at the right hand side. 

We now assume liquid helium to contain elementary excitations. For nonvanishing tem- 
peratures T, these excitations can be of thermal origin and thus give rise to the normal fluid, 
while at zero temperature they must be created by an external probe and yield a vanishing 
normal density p n . If T is above 1 K, the normal component is mainly a gas of rotons, being 
the phonons the dominant excitations at lower temperatures. Therefore, at any temperature 
the interaction of the vortex line with these elementary excitations, produces damped mo- 
tion of the vortex. As stated in the introduction, the main goal of this article is to construct 
a hamiltonian model which enables us to obtain this dissipative behavior. 

For this purpose, we shall consider a description of quantum dissipation similar to that 
recently presented in order to account for the irreversible evolution of solitons |Z3J], in which 
an effective Hamiltonian is constructed for the collective motion coupled to the residual 
excitations using the Collective Coordinate formalism |2"5|. This model exhibits unexpected 



features |26[ due to the fact that both the system and the reservoir have the same microscopic 
origin, which is just the case here discussed. 

In this spirit, and considering that within a superfluid in its ground state the vortex ex- 
hibits a soliton-like behavior, we propose a vortex-plus-reservoir Hamiltonian, that modifies 
expression (|2.1|) as follows 

H = ^ [p - qA(r) - A B] 2 + MSlv.y + H B (2.5) 

with B a vector function of operators that represents the elementary excitations of the 
superfluid and Hb is the Hamiltonian of these excitations. In Eq. ( |2.5| ) the interaction term 
Hint — — A B • v couples the reservoir and the vortex through the unperturbed velocity of 
the latter, being 



In the present approach, the hermitian operator B is associated to the creation of a 
density fluctuation in the liquid and could then be labelled by a transferred momentum q. 
Up to lowest order, one may have for each component of the vector B 

where O q (O q ) is the Feynman-Cohen operator that creates (destroys) a density fluctuation 
quantum, i.e., a phonon or a roton 

O q = P q -^E^/U- k - (2-8) 

being here N the number of atoms in the liquid. Furthermore, we realize that the term 
A 2 B 2 /2M appearing in (2.5) can be absorbed into the hamiltonian Hb, which is in charge 
of providing the equilibrium density vector of the reservoir. 

III. THE GENERALIZED MASTER EQUATION 



The Hamiltonian Q2.5|) is of the form system - plus - reservoir - plus - interaction |2T 



The standard reduction - projection procedure of nonequilibrium statistical mechanics [20 



enriched with the time convolutionless method developed by Chaturvedi and Shibata p2 
has already proven to be useful to derive a generalized master equation (GME), with time 
dependent coefficients, for the density operator a of a particle interacting with a heat reser- 



voir in the weak coupling - nonmarkovian limit [ 27, 28 1 . In this case, a is the density operator 



of the vortex and the generalized master equation reads |26 
da i A 2 ft 

— + -[H v , a] = --J J dr {[v x , [v x (-t), a]] + [v y , [v v (-t), (t]]}</>(t) 
A 2 ft 

-z-2 / dr {[v x , [v x (-t), cr] + ] + [v y , [v„(-t), a\ + ]} ip(r) (3.1) 



r Jo 

where [a, b] + denotes an anticommutator. In this expression, the time dependent functions 
4> and ip are the real and imaginary parts, respectively, of the correlation between heat bath 



operators fl29| , 



< B 3 (t) Bj >= 0(r) + iif)(r) (3.2) 

for j = x, y, assuming an isotropic reservoir. If the hermitian operator Bj is chosen according 
to Eqs. O, M, the function 

5 q (T) = q (r)+z^ q (r) (3.3) 

is just the Fourier transform of the dynamical structure factor S^q, uo) of helium II and is 
experimentally known for a wide range of transferred momenta |30|j3~i|| . 

Notice that the GME is a differential, rather than an integrodifferential, equation, since 
the unknown a under the integral sign is taken at time t; accordingly, it can be simplified if 
we define the following time dependent parameters, 

M A 2 r f 

- 7 (t) = -_ / dT^{T)smflr (3.4) 



2 ' w h Jo 
M , , A 2 r f 



M A r 

— fi(t) = — / dr mr) cosflr (3.5) 
2 n Jo 



and 



A 2 /■* 

Cm = -o / rfr0(r)cos^r. (3.6) 
n Jo 

The velocities appearing in Eq. ( |3.1| ) are those of the free vortex displayed in ( |2.6| ) and 
their detailed time dependence is extracted from Hamilton's equations corresponding to the 
Hamiltonian ( [2.1|) , namely 

v x{t) = [v x (0) — v s ] cosflt — v y (0) sinflt + v s 

v y (t) = [v x (Q) — v s ] sinfit + v y (0)cosflt. (3.7) 

In terms of these quantities and using Eqs. (3.4) to (3.6), we can write 

+ faff, ^} = -C(t) {[v x , [v x , a}} + [v v , [vy, a}}} (3.8) 
i M'yCt) rr r . . r r . 
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The effective Hamiltonian contains a renormalization to the vortex mass, induced by the 
coupling to the thermal reservoir, plus a drift contribution. Its expression is 

H eff = H V + (y 2 x + v 2 y ) + M v s w(t) v x — M v s j(t) v y (3.9) 

where u(t) = fi(t)\ Q=Q — fi(t). It is also worthwhile noticing that, being the system trans- 
lationally invariant on the (x,y) plane, terms in H e ff proportional to v x ,v y play no role in 
the dynamics. 

It is important to observe that the validity of the GME ( |3.8| ) is more general than the 
weak coupling approximation case. Indeed, if one expands the integral, time dependent 
collisional kernel of the master equation in powers of the coupling parameter A, as done, for 
instance, in Ref. |22[| , after a lengthy calculation one can realize that the form of the new 
GME is identical to ( |3.8| ), at least up to the fourth order in the expansion parameter, except 
for the fact that the coefficients (|3.4j) to (p.6|) become polynomials in A. 



Finally, we should also mention that the most common assumptions considered in many 
applications are that the reservoir is purely harmonic and/or that the interaction term is 
linear in its coordinates. If this is not the current case, the nonvanishing mean values < Bj > 
must be considered |26| and modify the effective Hamiltonian (|3.9|) . However, the equations 



of motion that we shall derive in the next section remain invariant, since these extra terms 
can be removed by a Galilean transformation. Moreover, it should be noticed as well that 
the correlation function < B^t) Bj > for k ^ j, which is in general a nonvanishing function, 
does not enter Eq. ( |3.1| ). 

IV. THE EQUATIONS OF IRREVERSIBLE MOTION 

We are now in a position to derive equations of motion for expectation values (a) of 
arbitrary observables a, which can be cast in the form 

'' + ~ < [a,H eff ] > = -C(t) (< [[a,v x ],v x ] > + < [[a,v v ],v v ] >) 



dt h 



+ £ (< [[a, v x ],v y } + > - < [[a, vy], v x } + >) . (4.1) 
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In order to derive equations of motion for the position and momentum components of the 
vortex we will restrict ourselves to the markovian limit; in other words, we consider that the 
correlation indicated by <f), ip is short lived, within the observational times. The parameters 
in Eqs. ( |3.4| ) and ( |3.5| ) become then time independent and after some algebra, elimination of 
the momentum permits us to write a unique complex differential equation for the expectation 
value of its velocity, that exhibits the effects of the coupling to the reservoir. This equation 



is 

72 



d 2 {z) n d(z) \ 

where the quantity that renormalizes the complex Magnus force is 

/3 = l + /i + z 7 . (4.3) 

with /i and 7 being the asymptotic values of ( |3.5|) and ( |3.4| ) respectively. 

We here realize that the reservoir constituted by the excitations of the superfluid provides 
both a dissipative and a conservative coupling, respectively measured by the parameters 7 
and /1. This is in agreement with the structure of the mutual friction force of the two fluid 
model JTJ-0. Moreover, keeping in mind that if the density fluctuations of the superfluid 
carry a definite momentum q, the heat bath correlation (y3 q y3 q (r)) is, to lowest order, just 
the Fourier transform of the dynamical structure factor, using standard relations of linear 
response theory |32|] one can readily show that 



where is the dynamical susceptibility or response function of the liquid (per unit 

length) at momentum q and energy HQ. In the most general situation where thermal ex- 
citations cover the whole momentum range, a summation over q should be applied on the 



right hand side of Eq. ( |4.4| ). A consequence of this result is that the temperature depen- 



dence of the drag coefficients is provided by the variation of x(q, fi) with T it is then 



worthwhile to keep in mind that the harmonic oscillator heat bath employed in previous 
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investigations of vortex coupling to excitations HTT| , |T3f predicts a temperature independent 
dissipation strength J29|. However, if the reservoir operators Bj are described by nonlinear 
functions of O}., O q rather than by the Feynman- Cohen operator ( |2.8| ), it is possible to show 
that the drag coefficients vanish at zero temperature |24|]2(J; in such a case, the coupling is 
uneffective and the vortex moves freely governed by the Hamiltonian ( |2.1| ), as expected. 

Equation ( t4.2|) can be straightforwardly integrated, giving a mean value of the complex 
velocity operator 



d(z) 



d{z) 



dt 



- v. 



t=0 



(4.5) 



Since Q 7 is always a positive quantity, this expression contains exponential decay of the 
initial conditions, and the limiting value of the vortex velocity is then the superfluid one v s . 

Let us now examine the situation as described by the phenomenological theory [l|-|3| 
where the drag force is written as 

dz \ 



Id = -(7o + ho) 



dt 



(4.6) 



where 70, 7 are the strenghts of the dissipative and conservative components and v n is 
the normal fluid velocity. It is also assumed that when equilibrium is reached, fo can be 
expressed in the form (a — ia') p s qh {v n — v s ). If one solves Newton's equation for a point 
particle with mass M moving under the Magnus and the drag force ( |4.6| ), one finds 



dz 
~dt 



dz 



dt 



- v s - (a' + i a) (v n - v s ) 



t=o 



o i{n-i iiM)t- lQ it/M 



+ (a +ia) (v n - v s ) + v s 



(4.7) 



where a and a' can be written in terms of 70 and 7g as 

To 



a = p s qh 



(p a qh-i ) +7o 



a 



/ 7o + 7o Hi -psQh) 



(4.8) 



(p s qh -i Q f + 7^ 

The inverse relationships giving 70, 7o in terms of a, a' can be found in Ref. 0. We see that 
in this case, the asymptotic velocity contains both the reactive and the resistive coefficients. 
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However, measurements of second sound attenuation in helium II at temperatures below 1.5 
K give values for a, a' around 10~ 2 , providing thus a negligible correction to the unperturbed 
velocity v s . 

On the other hand, it is important to notice that according to the general equation (4.1) 
derived in this work, the expectation value of the free velocity operator (2.6), whose complex 
counterpart reads 

v = JL +i 9. z (4.9) 
M 2 V ; 

satisfies the evolution law 



(v(t)) 



(«(0)> - «. I 1 - , 



e^ + vjl-^j. (4.10) 



with Hq = fiQ=o (cf. Eq. (3.5)). Equation (|4.10|) is remarkably close to the above expression 



( |4.7| ) for vanishing normal fluid velocity. Indeed, for a normal fluid at rest, Eq. Q4.10] ) is of 
the form (|4.7|), with coefficients 7o,7 (or a, a'), given by 

^- = a' + ia. (4.11) 
Elimination of /i and 7 gives the relationship 



a ~2 1 / 1, ~/\2 



a = -^r— ( 1 — ^TTi \m ( 4 - 12 J 



p s q h 70 

7o + (PsQ h - To) 
Ps q h (70 - p s qh) 

il + {p s qh-%f 

The interesting similitude between these relations and those in Eqs. (|4.8|) gives support 
to the conjecture that the present model embodies substantial aspects of the mechanism 
responsible of damped vortex motion in superfluid helium. The differences between the 
relationships characterizing the phenomenological model in Eqs. (4.8), and the present ones 
in Eqs. (4.12), are due to the fact that the structure of the drag force is not identical in both 
approaches. In fact, a close look at Eqs. (4.2), (4.3) and (4.6) shows that the Hamiltonian 
description gives rise to an extra component of the force, proportional to the relative two 
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fluid velocity v n — v s . This supplementary component is not removed by the assumption 
that the normal fluid lies at rest; however, it is also worthwhile to keep in mind that the 
assumption that the drag force is proportional to v - v„ applies under the hypothesis of 
vanishing vortex mass 0. If one is interested in getting rid of the extra force, a different 
model should be selected, so as to bring the two fluid dynamics into the picture. Such an 
improvement does not consist of a simple modification of the Hamiltonian (2.5); instead, 
a totally different formulation is required stemming from a Hamiltonian description of the 
two fluids to which a suitable copling is incorporated. This philosophy fits more specifically 
the spirit of macroscopic, fluiddynamical models and is thus beyond the scope of the present 
work. 

As a final remark, we wish to recall that every time dependent quantity here presented 
owes this dependence to the special model feature that makes room to a finite, although 
unknown, inertial coefficient of the vortex. This parameter rules the evolution since it 
appears in both the conservative and the decay time scale (cf. Eqs. (|4.5|), ( fL7| ) and ( |4 . 1 0| ) ) ; 
we then realize that as pointed out in Ref . |TD[ , experimental detection of the time dependent 
regime would thus provide a means of measuring the vortex inertia. It should be kept in mind 
that the present results concerning the dynamics cannot be extrapolated down to M = 0; 
in fact, the free Hamiltonian (2.1), the frequency (2.3), the vortex - reservoir Hamiltonian 
(2.5) and the velocity (2.6) become meaningless in such a case. However, the asymptotic 
velocity does not depend upon the mass, since its value causes the Magnus and the drift 
force to cancel each other in the absence of inertial effects. 

V. DISCUSSION AND SUMMARY 

Let us now examine further the characteristics of the model here presented and its 
relationship to the phenomenological description of dissipation. On the one hand, it is 
important to keep in mind that the vortex mass is assumed to vanish in the phenomenological 
two - fluid model of mutual friction, where the velocity arises from the balance between the 
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Magnus and the drag forces; it should be noticed that this regime is also the time - asymptotic 
form of a Newton - like equation of motion if the vortex mass is finite || . The precise value 
of the vortex inertial coefficient is thus not important in the limiting regime, although it 
influences the dynamics at finite times through the frequency Q (cf. Eqs. (2.3) and (|4.7|) ), 
which is the relevant parameter of the model. In this context, it is important to keep in 
mind that the coupling to the thermal excitations further renormalizes the vortex mass; in 
fact, inspection of the effective Hamiltonian (|3.9|) shows us that the kinetic energy has been 
changed into M(l + /i)v 2 /2. 

On the other hand, the phenomenological theories introduce a mutual friction whose drift 
and dissipative components are proportional to the relative two - fluid velocity v n — v s . The 
vortex velocity, either with respect to the superfluid or to the normal one, is determined by 
the force balance when inertial effects vanish; consequently, it depends upon the parameters 
of the drag. Instead, the present model should be regarded as a description in the reference 
frame of the normal fluid, i.e., both v s and d(z)/dt refer to the local velocity v n of the 
heat reservoir in the neighborhood of the vortex. The coefficients that measure the drag 
effects thus depend upon the strength of the coupling to the thermal excitations and upon 
their dynamical response; however the fluid dynamics of the elementary excitations is not 
explicitly contemplated. 

We believe that the model here presented covers most aspects of the description of 
dissipative dynamics of a vortex line in helium II and opens possibilities towards further 
improvements, among which, a definite one is the introduction of the motion of elementary 
excitations, to properly account for mutual friction in the sense of phenomenological theories. 
With respect to previous calculations of the drag coefficients carried, for example, in Refs. 
H||, our model, being quantal in nature, is not subjected to either the low - temperature 
limitations of a hydrodynamical description as pointed out in 0], or to uncertainties asso- 
ciated to a classical approach to the roton - vortex collisions [[J . It may be also mentioned 
that the model holds as well for vortex motion in liquid 3 He; quantum statistics only enters 
the characterization of the excitations making the heat reservoir, which would consist of the 
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zero sound phonons of the fermion liquid. No special differences with the present results 
would be expected in that case, except from the fact that the larger core size of vortices 
in 3 He could probably enlarge the inertia parameter M with a subsequent decrease in the 
oscillation frequency Q. 
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